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The unique property of Coulomb interaction in strict one-dimensional (ID) system is revealed 
that the Coulomb repulsion energy of paired electrons is divergent. As consequences, electrons in 
ID system can not doubly occupy the same spatial orbital and are completely localized. Numerical 
simulation by time dependent Hartree-Fock approximation shows this distinct property. The '0.7 
anomaly' in ID electron gas is fully explained by the property. Its possible contribution to strongly 
correlated systems is discussed. 



ID electron systems often demonstrate the strongly 
correlated characteristics, Coulomb interaction between 
electrons plays an important role but its essence still 
keeps elusive. For instance, in ID electron gases, there is 
a famous '0.7 anomaly'— phenomenon remaining unclear 
and controversial. It's regarded that the '0.7 anomaly' 
should be due to spin polarization and Coulomb interac- 
tion plays a key role^, but people do not understand why 
Coulomb interaction can lift the spin degeneracy. There 
is a latent need to illustrate the essence of ID Coulomb 
interaction and its affection on electron correlations. 

It is well known that density functional theory (DFT) 
loses its exactness for strongly correlated systems. In the- 
oretical field, Hubbard model^ is still widely employed 
for qualitatively explaining some strong correlation phe- 
nomenons. The key conception of Hubbard theory is the 
in-site Coulomb repulsion between paired electrons with 
opposite spins occupying the same spatial orbital. 
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A puzzling fact then rising up is that the Coulomb in- 
teractions have already been exactly and sufBciently in- 
cluded in the ab initio formalism of Hartree-Fock^ or 
DFT— method, considering the Coulomb repulsions not 
only between the paired electrons but also between all 
the other electrons in the other orbitals (or levels). In 
the closed-shell case, the Coulomb part of the Hamilto- 
nian in a single-particle equation if4 
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{V'jjJ = 1, 2. ..N/2} is the singular orbitals set, the 
Coulomb interactions from all the orbitals arc already in- 
cluded. If expanded by the atomic basis functions {(/)^} 
to generate the matrix representation, this operator is ex- 
pressed by the sum of a series of multi-center integrals^ 
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dridr2, (3) 



which is similar to the U in Eq. [U but all the other 
three and four-center integrals are there, means an evalu- 
ation of Coulomb repulsion more sufficient than Hubbard 
model. Therefore, it is very controversial to say that an 
ab initio computation is inclined to give an underesti- 
mated Coulomb interaction. On the other hand, when 
applying Hubbard model, the U can not be evaluated 
exactly but used only as an adjustable parameter. One 
possible reason for this controversial issue is that the ID 
Coulomb interaction has its distinct property which can 
not be embodied in conventional ab initio formalism. 

Some hidden presuppositions in the conventional ab 
initio formalism of Hartree-Fock, DFT or Hubbard model 
should be contemplated, one of them is that electrons 
doubly occupy orbitals. That presupposition is suspi- 
cious in ID electron systems, the property of Coulomb in- 
teraction determined by ID space dimensionality will be 
considered carefully in this article. Although literature^ 
proposes low dimensional Coulomb potential with a form 
different from 1/r in order to preserve the low dimen- 
sional Gauss's law, it is still assumed here that the po- 
tential keeps l/r form for ID systems, the ID electrons 
in the real physical world are just considered as confined 
three-dimensional (3D) ones. The basic idea in this arti- 
cle is that the property of the ID Coulomb integral in Eq. 
[T]or Eq. 13] is different from the 3D and two-dimensional 
(2D) cases. The ID space can have its unique topolog- 
ical property determined by its dimensionality. The pa- 
rameter U, or multi-center integrals in Eq. [3l should be 
deliberated according to the topological property. To do 
the two-electron (2e) integration in Eq. [1] it is common 
to decompose^ it into two steps, integration over ri and 
then integration over r2, that physically means average 
of Coulomb energy over electron 1 and then average over 
electron 2. Consider the first step. 
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The position at ri = r2 is usually regarded as a sin- 
gularity, especially in numerical simulations^. However, 
the singularity retains only for ID cases. For 2D and 
3D cases, it is actually not a singularity. We can employ 
polar (2D) and spherical (3D) coordinate systems to illus- 



2 



trate that point. Shift the origin to r2, correspondingly 
replace |ri — by r and ^(ri) becomes (j){r,uj), uj is 9 
for 2D and {6, (/>) for 3D, then the integral is rewritten as 

/ = J f^^n^r--'dr. (5) 

Here n is the dimensionality. There is no singularity 
for 2D and 3D cases (n > 1) because the denominator 
r can be compensated by the integration volume ele- 
ment r"'~^dr. This property determines whether the 2e 
Coulomb integral of Eq. |3]is infinite or not. There are 
already several works^"— about the 2e integrals for 2D 
systems. Using Slater- typoji^iii or Gaussian- typc^ ba- 
sis functions, analytical results for 2D Coulomb integrals 
can be obtained and they are finite, as expected. Numeri- 
cal simulationi^ by time-dependent Hartree-Fock method 
also shows the numerical Coulomb integral between two 
partially superposed 2D electron waves is limited. As 
regards in 3D cases, the formulas of 2e integrals using 
the Gaussian-type basis is well-built^— and widely used. 
From their deduction of these formulas we can see that 
the compensation by the r"~^dr plays a key role to make 
the singularity disappear and then integrals finite. The 
point is, that finitcncss could not hold for ID system 
and therefore, the ID 2e integrals is divergent if there is 
any superposition of wave functions. ID electron's be- 
havior should be quite distinct and we need consider it 
from a generalized viewpoint. To my knowledge, most of 
ID system modeling works only escape that singularity 
by using soft-Coulomb potentia l^^'^^ , i.e. to cut off the 
Coulomb potential, or replacing it by i5 potentiali^. This 
simplification by a cutoff could not give results with real 
physical meaning, it's shown below that different cutoff 
value can remarkably change the electron's behavior. 

Actually, in some case we can have the analytical form 
of ID 2e Coulomb integral and easily show it is diver- 
gent. Consider the ground state of a ID parabolic po- 
tential well, the spatial orbital is a Gaussian wave packet 
■0 = \fCe~'^ suppose it is occupied by two electrons 
with up and down spins. The two coordinates of xi and 
X2 actually form an infinite plane, a transform to the 
polar coordinate system, xi ~ rcosO^ X2 = rsinO and 
dxidx2 rdrdO, then can be applied. The Coulomb 
integral becomes 

^^^dx,dx2 = / , . 

\xi—X2\ Jo \cosa — stny\ 

(6) 

This integration is divergent while 9 —> tt/A and 9 — > 
57r/4, where \xi — a;2| — >■ 0. If the Hubbard model is 
applied here, that means the parameter U in Eq. [T] is 
infinite. An infinite Coulomb energy is non-physical and 
does not exist. So, it is very naturally to conceive that in 
a completely strict ID system, electrons cannot occupy 
an orbital in pair. Moreover, no superposition of wave 
functions of different electrons is permitted because that 
will result in singularity in the Coulomb integrals and 
divergent energy, hence ID electrons arc always localized. 



To show this distinct property's effect on the behavior 
of ID electrons, a numerical simulation is done to two 
electrons in a ID infinite potential well. Time-dependent 
Hartree-Fock approximation^^ is employed. The two 
electrons are assigned opposite spins in order to neglect 
the exchange interaction to simplify computation. Then 
two time-dependent single-particle equations are built as 

ih = -— ViVi(a;i)+- / dx2- :^i{xi) 

ot 2m 47reo J \xi — X2\ 

(7) 

.^dip2{x2) -^^^2, / ^ /"j i^H^i) , , ^ 

ih = ^T—^ 2i'i(x2)+- / dxi- rV'2 X2 

at 2m 4:7160 J \xi — X2\ 

(8) 

The last item in each equation is the Coulomb repulsion 
potential from each other. The equations are discretized 
by the finite difference time dependent (FDTD)ii ap- 
proach, the Coulomb integrals are calculated by numer- 
ical integration method. We face the problem of how 
to deal with the singularity in the numerical integra- 
tion. The ID space is discretized into many grid points 
{Xi}, the singular points are those where xi = X2, i.e. 
I^^i — 3^2 1 = 0- One choice is to ignore the singular points 
while doing numerical integrations but make the spatial 
difference AX = \Xi-Xi±i\ very small. When AX -> 0, 
the integral will approach to its real value. This method 
is tested effective for 2D and 3D numerical integrations 
in which a converged finite integral value very close to 
the real value can be obtained even when AX (Ay, AZ) 
is not necessary to be super small. In the ID case, the 
integral is divergent when AX — > 0, but the speed is not 
quite fast, that means a very fine AX is needed and it 
hugely enhances the time costs. We can make an other 
choice to avoid expensive time costs resulted from too 
fine AX. The singular points are included, however, in- 
stead of set \xi — X2\ as 0, a finite but very small value is 
to make a substitution of it. In this simulation, the value 
is 10~^A. This is actually a cutoff to the Coulomb repul- 
sion potential. For numerical simulation it is reasonable, 
since we only want to sec the effect of a large enough 
Coulomb repulsion. But keep in mind that only a small 
enough cutoff can give you the real physics picture. 

The initial wave functions are set as two Gaussian wave 
packets defined as 

V'i,2(a;, 0) = e-"(^±^o)=+*fei,2a. (9) 

In my simulation, a ~ lOOA ^, .tq = 0.5A, fci = 
lOOA . k2 ~ — lOOA . they arc moving to each other 
to generate a collision. The value of k indicates a high 
kinetic energy, about 41keV, ensuring the collision very 
strong. Fig. 1 shows the collision process and the evolu- 
tion of the wave functions. From t = 32000AT, the first 
collision starts, it's finished at about t ~ 70000AT. Dur- 
ing the collision process, no superposition takes place. 
Electrons totally block each other and no tunneling 
through is permitted. As time passes, the wave pack- 
ets just collide again and again and finally disperse into 
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FIG. 1. (Color online) The collision process and evolution 
of the two electron wave packets in an infinite potential well 
simulated by FDTD method. The width of the well is 2 A. 
Space difference is AX = O.OOf A; temporal step is AT — 
10~* femtosecond, (a) is the initial state; from (b) to (e) 
show the first collision process; (f) is the state after a long 
time of evolution, showing the dispersion and localization of 
the electron wave packets. 



many chaotic peaks in the well. The dispersed waves 
keep localized in their original sides. They do not form 
a pair to share the same orbital as they are supposed 
to do. One may argue that in Fig. 1(f), some tiny su- 
perpositions do appear here. That is because we use a 
finite cutoff to approximate the infinite Coulomb repul- 
sion and they arc just numerical errors. Actually, if the 
cutoff value is not so tiny enough, wave packets will pass 
through each other. That shows that an arbitrary cut- 
off of the Coulomb potential may not give you the real 
physics picture of ID system. One can imagine that for 
real Coulomb repulsion, i.e. an infinite small cutoff, this 
superposition should disappear completely, that is right 
the real physics in an ideal ID system. 

The conclusion that no superposition of wave functions 



is permitted in an ideal ID system is extensively mean- 
ingful. Concisely, two ID electrons with opposite spins 
can not occupy the same spatial orbital as they do in 2D 
and 3D space. More strictly, ID electrons are totally lo- 
calized. Neither electrons pass through each other nor do 
they overlap with each other at all. Consequently, the ar- 
rangement order of electrons is determined and electrons 
arc distinguishable since there is no any superposition (or 
overlap). In other words, Pauli repulsion disappears in 
ID system. Consider the exchange operator Ka 
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If ipa and ipi have no superposition, the integration is 
zero, then the exchange energy is also zero. For distin- 
guishable particles, Fermi statistics is not necessary and 
non-applicable. In view of the facts, we can draw a fur- 
ther conclusion that no band theory for ID electrons is 
applicable. Many textbooks of solid physics start from 
the ID periodic system to introduce the picture of band 
theory, but that is just not appropriate. Periodic wave 
functions of many electrons distributed in full ID space 
must induce superpositions, which consequently result in 
infinite Coulomb repulsion energy. ID electrons arc al- 
ways strongly correlated, free particle approximation is 
not applicaple. One might argue that for some particu- 
lar quasi-lD systems such as carbon nanotube or atom 
chains, the band theory works very well. It is because 
these particular systems are not completely strict ID but 
quasi-lD, their electron wave functions are 3D. A quasi- 
lD system is actually 3D one in which only the motion of 
electrons along one special direction is more prominent 
than the other two directions. Only when electrons are 
strongly constrained in one dimension, the ID property 
can take on and its perturbation could not be ignored. 

Completely strict ID system is the extreme limitation 
of a constrained 3D system. To some extent, a quasi-lD, 
i.e. constrained 3D, system is physically more realistic 
than a pure ID one. Olszewskii^ investigated such a 
system, interacting electron gases in a long but narrow 
cylindrical tube. When the radius of tube shrinks to zero, 
it degenerates to pure ID. Exact analytical solutionsiS. of 
free electron in 3D cylindrical box 

= CNi.stn{^)MXi^)e'''^ (11) 

are employed as the basis functions to perform conven- 
tional Hartrec-Fock calculation. Here, L is the length and 
i? is the radius of the tube. The radial part Ji{Xij^) is the 
Bessel function, the angular part is e*'"^. Notice that the 
longitudinal part sin{^^^^) is some sort of plane wave ba- 
sis function distributed in the whole z space, then the su- 
perposition of electron waves is supposed to be there. Us- 
ing these basis functions, both Olszewskii^ and Robertai^ 
show that the Coulomb interaction integrals go to diver- 
gent while taking the limit i? — > 0. This result is consis- 
tent with my discussion above. More interesting results 
are the total energies of the ground states of the gases 



4 



calculated by Hartrce-Fock method. As i? — >■ 0, the dou- 
bly occupied case has ground state energy much higher 
than the singularly occupied case, and both of them are 
divergent while i? — > 0. These results are not surprising. 
The pairing and superposition of electron waves in a pure 
ID system must result in infinite Coulomb energy. Infi- 
nite energy is non-physical, then we can imagine when 
the i? — )■ 0, i.e. the constraint becomes strong, the real 
physics picture should be that the electrons must depair 
and localize themselves to lower the energy. If one force- 
fully applies conventional Hartree-Fock formalism to an 
ideal ID system, only some non-physical results could 
be obtained because it forces electrons to be paired and 
delocalized. 

In ID electron gases, there is a long standing puzzle 
called as '0.7 anomaly'— which is often considered in- 
duced by strong correlation. This phenomenon can be 
easily explained by my conclusions for ID systems. In 
the experiment by Thomas et alX, a gate voltage Vg is 
applied to 2D electron gases (2DEG) to create ID con- 
strictions and control the structure width. When increas- 
ing the strength of Vg, the constriction becomes strong 
and the width of 2DEG becomes narrow, then finally 
the 2DEG is transformed to ID gas at a large enough 
Vg. During this procedure, the conductance channels are 
closed one by one and a series of plateaus generate. Each 
plateau corresponds to a quantized conductance le^ jh. 
The 2 is due to the doubly occupying on one channel. 
After the last plateau where the constriction is strong 
enough to squeeze the 2DEG to ID, the anomaly value 
0.7(2e^//i) appears. This is not surprising either. As dis- 
cussed above, when the constriction is very strong, elec- 
trons should access the ID limitation and take on their 
ID behaviors and depair themselves. A totally depairing 
makes each channel have only one electron, the quan- 
tized conductance should be 0.5(2e^//i). If the system 
can not totally access the ideal ID limitation, a partially 
depairing will happen and result in a conductance be- 
tween 0.5(2e2//i) and \{2e^ jh). Denote it as ^{2e^ jh). 
The coefficient 7 depends on the strength of constriction 
and the energy of conducting electrons. The value '0.7' 
is only one special case. Reilly et al^ measured long 
and clean ID wires. Very clearly the plateaus appear at 
0.5(2e^//i) while increasing the length. That's because 
long enough constriction makes the ID effect remark- 
ably prominent and the electrons almost totally depair. 
Thomas et al^ suspect there is spontaneous lifting of spin 
degeneracy in the ID constriction due to electron inter- 



actions. That guess is correct, however, the lifting of de- 
generacy should not be ascribed to magnetic mechanism 
but to the ID Coulomb mechanism. This '0.7 anomaly' 
is a very good evidence to support my conclusions for ID 
Coulomb interaction. 

Up to now, theoretical researches for strongly corre- 
lated systems, even for many quasi- ID problems, arc all 
based on these presuppositions that electrons occupy spa- 
tial orbitals in pair, they obey Fermi statistics, and un- 
questionably the band theory is applicable and applied. 
The idea that the spatial orbital be occupied by two elec- 
trons originated from early quantum theory motivated by 
experiments of that time and then is widely adopted in 
modern quantum chemistry and physics theories. The 
spatial parts of electron waves can be totally or par- 
tially superposed in 2D and 3D spaces. As the costs, the 
Coulomb repulsion energy is limitedly enhanced. How- 
ever, inside some special systems, if there are very pure 
ID structures or strong ID constraints resulted from 
pressure, electrical or magnetic forces, these presuppo- 
sitions are questionable. The property of ID Coulomb 
interaction forbids any superposition of electron waves. 
Very possibly it is the ID Coulomb interactions that 
make strong perturbations to these systems and make 
their behaviors so strange. Therefore, for these special 
systems, the formalism of conventional DFT based on 
3D basis functions and 3D Coulomb integrals calculation 
in Eq. [3]cannot get correct results. Under strong ID con- 
straints, the electron waves arc distorted, the Coulomb 
interaction is neither pure 3D nor pure ID. Wc need 
develop new methods to exactly calculate these tricky 
Coulomb interactions. Perhaps this is also one reason 
why it can give us results more reasonable when we use 
it only as an adjustable parameter in Hubbard model. 

In summary, the property of the ID Coulomb interac- 
tion is revealed. Electrons are completely localized in ID 
systems, they couldn't doubly occupy one orbital, and 
the Pauli repulsion does not exist. Surprisingly these 
properties have been ignored all along. Many theories 
and models for ID electron systems just naturally adopt 
the assumption of doubly occupying and use the Slater 
determinant form of wave function to include the ex- 
change effects. This is the first time to definitely con- 
clude that electrons in their ID limitation will depair 
and localize themselves. 
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